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Abstract 

We study a large family of competing spatial growth models. In these the vertices in 1 d 
can take on three possible states {0,1,2}. Vertices in states 1 and 2 remain in their states 
forever, while vertices in state which are adjacent to a vertex in state 1 (or state 2) can 
switch to state 1 (or state 2). We think of the vertices in states 1 and 2 as infected with 
one of two infections while the vertices in state are considered uninfected. In this way 
these models are variants of the Richardson model. We start the models with a single vertex 
in state 1 and a single vertex is in state 2. We show that with positive probability state 1 
reaches an infinite number of vertices and state 2 also reaches an infinite number of vertices. 
This extends results and proves a conjecture of Haggstrom and Pemantlc 5 . The key tool 
is applying the ergodic theorem to stationary first passage percolation. 

1 First Passage Percolation 

In this paper we study a class of competing spatial growth models by first studying stationary 
first passage percolation and then applying our results to the spatial growth models. In first 
passage percolation every edge in a graph is assigned a non-negative number. This is interpreted 
as the time it takes to move across the edge. This model was introduced by Hammersley and 
Welsh 0. See for an overview of first passage percolation. 

Let [i be a stationary measure on [0, oo)^^S es ( zd ) an d let w be a realization of \i. For any x 
and y we define the passage time from x to y, r(x,y), by 

r{x,y) = inf y~]uj(vi,v i+1 ) 

where the sum is taken over all of the edges in the path and the inf is taken over all paths 
connecting x to y. 

The most basic result from first passage percolation is the shape theorem. We let = 
(0, . . . , 0) and 1 = (1, 0, . . . , 0). Define 

S(t) = {x : t(0,x) < t} 



1 



and 



S(t) = S(t)+ 



The shape theorem says that there is a nonempty set S such that — ^ converges to S a.s. 

Theorem 1. ^ Let \x be stationary and ergodic, with the distribution on any edge have finite 
d + e moment for e > 0. There exists a set S which is nonempty, convex, and symmetric about 
reflection through the origin such that for every e > there exists a T such that for all t > T 



This theorem is a consequence of Kingman's subadditive ergodic theorem. It is the only 
property of first passage percolation that we need. In general little is known about the shape of 
S other than it is convex and symmetric. Cox and Durrett have shown that there are nontrivial 
product measures such that the boundary of S contains a flat piece |2|. However for any compact 
nonempty convex set S there exist a stationary measure [i such that the shape for /i is S [I]. 

Another widely studied aspect of first passage percolation are geodesies. A geodesic is a 
path G = {vq ,«!,...} such that 



for any m < n. We let G^{x,y) = G(x,y) be the union of all geodesies connecting x and y. 
Define 



We refer to this as the tree of infection of x. We define K(T(x)) to be the number of 
topological ends in T(x). This is also the number of infinite self avoiding paths in T(x) that 
start at x. 

Newman has conjectured that for a large class of \x, |if(r(0))| = oo a.s. [S] Haggstrom 
and Pemantle proved that if d = 2 and [i is the i.i.d with exponential distribution then with 
positive probability |ET(r(0))| > 1. Newman has proved that if /U is i.i.d. and S has certain 
properties then |iT(r(0))[ = oo a.s. jH] Although these conditions are plausible there are no 
known measures \i with S that satisfy these conditions. 

Now we will introduce some more notation which will let us list the conditions that we place 
on fj,. We say that the configuration u> has unique geodesies if for all x, y € 7L d there exists a 
unique geodesic from x to y. If there exists a unique geodesic between x and y we denote it by 
G(x,y). The configuration co has unique passage times for all x and y ^ z 




n-l 




T(x) = U yeZ d{e : e G Edges(Z d ) and e G G(x,y)}. 



i~(x,y) / t(x,z). 
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For any u) we let ^ ' be the conditional distribution of \i on the edge (0, 1) given that uj'(v, w) = 
co(v,w) for all edges except (0, 1). We say that \i has finite energy if for any set A C R such 
that /x{u;(0, 1) £ A} > and almost every u, /i^' 1 ' {lo' : u/(0, 1) G A} > 0. 

As /x is a stationary measure we can study its ergodic theoretical properties. For any v € T, d 
define the shift map T v : [0, oc) Ed g es ( z<i ) [0, oo) Ed g es ( z<i ) by 

T>)(j)=u,(j + <;) 

for all j € Edges(Z^). The measure fi is totally ergodic if for all v £ Z d the action (/x, T 1 ') is 
ergodic. 

Now we are ready to define the class of measures that we will work with. We say that \i is 
good if 

1. n is totally ergodic, 

2. fj, has all the symmetries of Z d , 

3. the distribution of \i on any edge has finite d + e moment for some e > 

4. /x has finite energy 

5. /ui 0,1 ^ is an absolutely continuous measure with support [0, oo) a.s., and 

6. /i produces a shape S which is bounded. 

Note that conditions [21 IU and [5] imply that fj, has unique geodesies and unique passage times. 
These conditions were chosen to make the arguments as easy as possible and could be made 
more general. All that is essential for the argument to show that there are at least two disjoint 
infinite geodesies is that ^ is totally ergodic and that Lemma ^ an d Corollary ^ below are 
satisfied. The conditions El \M an d IS) are used to show that coexistence occurs with positive 
probability. Throughout the rest of the paper we will assume that \i is good. Unfortunately 
there is no general necessary and sufficient condition to determine when the shape S is bounded. 
See [3J for examples. 

2 Spatial Growth Models 

Now we explain the relationship between first passage percolation and our competing growth 
models. For any uj € [0, oo) Ed § es ^ d ^ with unique passage times and any x ^ y £ Z, d we can 
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project it to Co x>y € {0, 1, 2} z by 





= < 



2 if t(x, z) < t and t(x, z) < r(y, z); 
1 if r(y, z) < t and t(x, z) > r(y, z); 
else. 



If fi has unique passage times a.s. then /x projects onto a measure on I {0, 1, 2} 



[0,oo) 



. It is clear 



that the models start with a single vertex in state 1 and a single vertex is in state 2. Vertices 
in states 1 and 2 remain in their states forever, while vertices in state which are adjacent to 
a vertex in state 1 (or state 2) can switch to state 1 (or state 2). We think of the vertices in 
states 1 and 2 as infected with one of two infections while the vertices in state are considered 
uninfected. In this way these models are variants of the Richardson model. 

As each z € Z rf eventually changes to state 1 or 2 and then stays in that state for the rest 
of time, we can speak of the limiting configuration. There are two possible outcomes. The first 
is coexistence or mutual unbounded growth. If this occurs then the limiting configuration has 
infinitely many z in state 1 and infinitely many z in state 2. The other outcome is domination. 
If this happens then in the limiting configuration there are only finitely many vertices in that 
state and all but finitely many vertices are in the other state. 

For many measures \x (for example if \x is i.i.d. with nontrivial marginals) then it is easy to 
prove that domination occurs with positive probability. But it is much more difficult to show 
that coexistence occurs with positive probability. More precisely we define C(x,y) to be the 
event that 



We refer to this event as coexistence or mutual unbounded growth. Our main result is 
that with positive probability coexistence occurs. 

Theorem 2. If fi is good then 



This proves a conjecture of Haggstrom and Pemantle 0. They proved this theorem in the 
case that d = 2 and fi is i.i.d. with exponential distribution. Garet and Marchand have given 
a different proof of Theorem |^1 1 31 . Their method follows more closely the approach taken by 
Haggstrom and Pemantle. 



\{z : lim lo x Jz) = 1}| = \{z : lim u xy (z) = 2}| 

i— s-oo t — >oo 



= 00 « 



P(C(0,1))>0. 



4 



3 Outline 



In this section we outline the proof of our main result. For any x,y £ Z rf and infinite geodesic 
G = (vq, vi, V2, ■ ■ ■ ) we can define 

b g( x i V) = b g{x, y) = lim t(x, v n ) - r(y, v n ). 

n— too 

To see the limit exists first note that 

B G (x,y) = lim r{x,v n ) - r(y,v n ) 

= lim r(x, v n ) - t(v , v n ) + t{vq,v u ) - r(y, v n ) 

n— >oo 

= lim (t(x, v n ) - t(v , v n )) + lim (t(v , v n ) - r(y, v n )). 

n—>oo n— >oo 

As G is a geodesic the two sequences in the right hand side of the last line are bounded and 
monotonic so they converge. Thus Bc(x,y) is well defined. If for a given u and all x,y € 1j d 
the function Ba(x,y) is independent of the choice of infinite geodesic G then we can define the 
Busemann function 

B(x,y) = B u (x,y) = B G (x,y). 

The main step in our proof is LemmaU which states that the probability that {B(x,y)} x yG %d 
is well defined is 0. 

We will work by contradiction to prove Lemma |IJ In Lemmas [5] and |21 we assume that 
{B(x, y)} Xj y € z d 1S weu defined a.s. and then apply the ergodic theorem to {B(x, y)} Xi y^zd- Then 
in Lemma 0] we show that the conclusions of Lemma 01 generate a contradiction with the shape 
theorem. Thus with positive probability there are vertices x and y and distinct geodesies Go = 
Gq{uS) and G\ = G\{u) such that 

Ba (x,y) / B Gl (x,y). 

From this point a short argument allows us to conclude that coexistence is possible with positive 
probability. 

4 Proof 

The heart of the proof is applying the ergodic theorem to the Busemann function. This is done 
in Lemmas |^1 and 01 We start by showing that the symmetry of \i implies that the expected 
value of the Busemann function is 0. 
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Lemma 1. If {B(x,y)} x y£ id is well defined a.s. then for all v € Z rf 

B(B(0,v)) = 0. 

Proof. By symmetry of /i we have that E(i3(0, 1)) = E(£?(l, 0)). Combining this with the fact 
that B(0, 1) + B(l, 0) = proves the lemma. □ 

Now we apply the ergodic theorem to B(0,v). 

Lemma 2. If {B{x,y)} x y&Ii d is well defined a.s. then for all v € 1* d and e > there exists M 
such that 

P(\B(0, mv)\ < em for all m > M) > 1 — e. 
Proof. First rewrite B(0,mv) as follows. 

B(0,mv) = B(0, v) + B(v, 2v) H h B((m - l)v, mv) 

B(0,mv) = B UJ {0,v)+B TV ^\0,v) + --- + B T(m ~ 1)v ^(0,v) 

m—l 

B(0,mv) = J2B Tiv ^(0,v) (1) 

As [i is good it is totally ergodic and the action (T v ,fx) is ergodic. Thus by line (^Q) and Lemma 
^the claim is a consequence of the ergodic theorem. □ 

We now strengthen this lemma by using the following corollary of shape theorem. For x € 1j d 
we let \x\ = \x%\ + \x2\ + • • • + \xd\- 

Corollary 1. There exist < k\ < k<i < oo such that for every e > there exists an N such 
that 

( t(0 x) \ 
P \ k\ < -. ' < fc 2 /or a// x such that \x\ > N j > 1 - e. 

Proof. The existence of k<i is due to the fact that the set S (from Theorem ^) is nonempty. The 
existence of k\ follows because one of the requirements of \i being good is that S is bounded. □ 

Lemma 3. If {B{x,y)} x y& id is well defined a.s. then for any e > there exists N such that if 
n> N then 

/B(0,x) \ 
P I — — - — < e for all x such that \x\ = n ) > 1 — e. 



Proof. Given e > pick vectors vi,v%, . . . ,Vj such that \v\\ = \v2\ = ■■■ = \vj\ and for all x 
sufficiently large there exists i S {1, 2, . . . , j} and m € N such that 

\x — mvA < e\x\ and mlvA < \x\. 
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For all x and y we have that 

B(0,x) = B(0,y)+B(y,x). 

This implies that for any x and y 

B{0,x)<B(0,y)+T(y,x). 

For any n let m be the largest integer such that m\vi\ < n. (This is independent of i.) Thus if 
there exists x with \x\ = n and — tj-^ > e then there exists i such that either 

1. B(0,mvi) > en/2 = e\x\/2, or 

2. \x — mvi\ < e\x\/2k2 and t(x, mvi) > e\x\/2. 

(The constant ki is from Corollary ^) 
By Lemma HI there exists M such that 

P(there exists m > M and i £ {1, 2, . . . , j} such that B(0,mvi) > 2em\vi\/3 > en/2) < 2e/3. 

Thus the probability of the first event is less than 2e/3 if n is sufficiently large. 
By Corollary ^ there exists L such that for any I > L 

P(there exists z with \z\ < I and r(0, z) > hit) < e/3j. 

Applying this with each mv{ in place of and en/2/c2 in place of I we get that the probability of 
the second event is less than e/3 if n is sufficiently large. Thus for any e > we get N so that 
if n > N we get that 

P (there exists x such that \x\ = n and — ^— - — - > e ) < e 

V M / 

which proves the lemma. □ 

Next we show that this generates a contradiction with the shape theorem. 
Lemma 4. P({B(x, y)} xy ^z d ^ s we ^ defined) = 0. 

Proof. We work by contradiction. Suppose that with positive probability {B(x, y)} x ^ y& i<i is well 
defined. The Busemann function being well defined is a shift invariant event which, by the 
ergodicity of fx, implies that {B(x,y)} x yeZ d is well defined a.s. and the conclusions of LemmaEl 
apply. Pick e < i min(A;i, 1), where k% comes from Corollary ^ By the choice of e and Corollary 
^we have that there exists N such that for all n > N 

/rfO x) \ 2 

P \ \ ' > 2e for all x such that \x\=n) > -. (2) 
V \x\ / 3 
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By Lemma El there exists n > N such that 

B(0,x) \ 2 

— — - — < e for all x such that \x\ = n I > — . (3) 
\x\ J 3 

However there exists at least one infinite geodesic G = (0, vi, V2, ■ ■ ■ ) which begins at 0. (The 
choice of G is immaterial.) For all n there exists k such that \vk\ = n. For any k we have that 
B(0,Vk) = r(0, Vk). This shows that lines ((2J) and © cannot both be true. Thus the lemma is 
proven. □ 

Note that the lack of a well defined Busemann function implies that there exists at least two 
disjoint infinite geodesies. Now we show that the lack of a well defined Busemann function also 
implies coexistence has positive probability. Coexistence is implied if there exist two infinite 
geodesies Go = (vo,v\,V2, • • • ) and G\ = (u>o, w\, W2, • • • ) such that 

B Go (0,l)<0< B Gl (0,l). 

We show coexistence is possible by showing that we have two such geodesies with positive 
probability. 

Proof of Theorem [2} By Lemma 0] we get an event A of positive probability and x,y £ TL d 
such that for all u> € A we have two geodesies Go = Go(cj) = (uq, v\, V2, ■ ■ ■ ) and G\ = G\{uj) = 
(w ,wi,W2, • • • ) with 

Bg q (x,V) < B Gl (x,y). 

(If there is more than one pair of geodesies which satisfy this equation we can choose Go and 
Gi in any measurable manner.) It causes no loss of generality to assume that \x — y\ = 1. Thus 
by the symmetry of \i we can assume that x = and y = 1. As B Go (0, 1) and Bq x (0, 1) do not 
depend on any finite number of edges in the geodesies, it causes no loss of generality to assume 
that 0, 1 are not endpoints of any of the edges in Go or G\. By restricting to a smaller event 
A C A of positive probability we get a nonrandom r > such that for all uj £ A 

B Go (0,l)<r<B Gl (0,l). (4) 

By the symmetry of [/, we can assume r > 0. From the definition of B Gl (0,l) we get that 
B Gl (0,l) <r(0,l). 

Now we form a new event A'. Given u> £ A define u 1 by 



oj'(v, w) 



uj(v, w) + r, if 1 € {v, w}; 
uj(v,w), else. 
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The event A' consists of all uj' that can be formed from in this way from some to E A. By 
conditions |2 and El of the definition of jj, being good, the event A' also has positive measure. 
We let t' indicate the passage times in ui' and r indicate the passage times in w. It is easy to 
check that for any z / 1 

r (1, z) = t(1, z) + r. 
Also if 1 is not an endpoint of any of the edges in the geodesic G w (0, z) then 

t'(0,z) = t(0,z). 

As Bg (0, 1) < Bg±(0, 1) < t(0, 1) we have that for all large n the vertex 1 is not an endpoint 
of any of the edges in the geodesic G w (0, v n ). Thus t'(0, v n ) = r(0, v n ) for all large n. Also note 
that since neither or 1 is an endpoint of any of the edges Gq{uj) or Gi(ui) we have that Gq(uj) 
and Gi(uj) are both geodesies for uj' . 
Thus for any uj' € A' we have that 

B^ oH (0,l) = lim o r , (0,v n )-r'(l,v n ) 

= lim t(0, v n ) - (r(l, v n ) + r) 

n— »oo 

< 0. 

The last step follows from line (@J). We also get that 

B# lM (0,l) = lim T'(0,w n )-T'(l,w n ) 

> lim r(0, iu n ) - (r(l, w n ) + r) 

> BS lM (0,l)-r 

> 0. 

The last step follows from line @. Thus we have coexistence for all w' £ i'. 
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